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ABSTRACT: Using an analytic model we show that the predictions of the cosmic string 
model for the peculiar velocities and the microwave background (MBR) anisotropy depend 
on similar combinations of string evolution parameters. Normalizing from the COBE 
detection of MBR anisotropy, and for certain reasonable values of string network evolution 
parameters, we find that the magnitude of predicted velocity flows is in good agreement 
with observations on small scales but is inconsistent with observations on large scales (> 
50/i -1 Mpc). The Cosmic Mach Number obtained from the cosmic string scenario is found 
to depend on a single network evolution parameter and is consistent with observations on 
scales 5fr _1 Mpc to 20/i _1 Mpc. 



The cosmic string model for large-scale structure formation with hot dark matter 
promises to be successful in many ways (Vachaspati T. & Vilenkin A. 1991; Vollick D. 
N. 1992; Albrecht A. & Stebbins A. 1993). The scenario might automatically account for 
large-scale filaments and sheets (Vachaspati T. 1986; Stebbins A. et al. 1987; Vachaspati 
T. & Vilenkin A. 1991; Perivolaropoulos L., Brandenberger R. & Stebbins A. 1990; Hara 
T. & Miyoshi S. 1990), galaxy formation at epochs of z ~ 2 — 3 and galactic magnetic fields 
(Vachaspati T. & Vilenkin A. 1991; Vachaspati T. 1992b). It can also provide (Vachaspati 
T. 1992a; Brandenberger R. et al. 1987) large-scale peculiar velocities (Dressier A. et. al. 
1987; Burstein D. et. al. 1983; Peebles P. J. E. & Silk J. 1990) and is consistent (Bennett 
D. P., Stebbins A. & Bouchet F. R. 1992; Perivolaropoulos L. 1993a,b) with the measured 
anisotropy of the microwave background radiation (Smoot G. et. al. 1992). 

Even though the cosmic string model involves a single free parameter (the mass per 
unit length of the string), its predictions depend also on our understanding of the evolution 
of the string network. The main features of this evolution can be encoded into a set of 
parameters which have been constrained by recent detailed numerical simulations (Bennett 
D. P. & Bouchet F. R. 1988; Allen B. & Shellard E. P. S. 1990) (see also papers by Albrecht 
& Turok in (Gibbons G. W., Hawking S. W. & Vachaspati T. 1990)). Therefore, in 
order to test the cosmic string model, it is important to obtain expressions of the model's 
observational predictions in terms of these parameters and then use observations to fix 
their values. The consistency of the model can then be tested by comparing the fixed 
values with numerical simulations of string evolution and with constraints from different 
observations. 

In this letter we show that the predictions of the peculiar velocities and MBR anisotropies 
resulting from cosmic strings depend on roughly the same combination of undetermined 
parameters which we call a. The prediction of peculiar velocities also depends on pa- 
rameters denoted £ and ^5/4, denoting the curvature radii of strings and the Newtonian 
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interaction of wiggly strings with matter. On the other hand, the microwave background 
predictions are independent of / 5 / 4 (since there is no Newtonian interaction of strings with 
photons) . By combining the two predictions and comparing with observations, we can im- 
pose constraints on these parameters, thus providing an interesting test for the cosmic 
string model. One outcome is that the string scenario cannot account for the observed 
peculiar velocities on all the different scales. With the most suitable choice of parame- 
ters, there is good agreement with observations on small scales but inconsistency with the 
observations on large scales. This potential problem of cosmic strings however, is also en- 
countered in several other models (Brandenberger R. et. al. 1987) including the standard 
CDM model. 

It is not difficult to guess that the peculiar velocities and the microwave background 
anisotropy should depend on roughly the same combination of network parameters. The 
peculiar velocities arise because with every string sweeping across the horizon, any fixed 
volume of matter experiences an impulse towards the wake of the string. At the same time, 
the anisotropy arises because the string gives an impulse to the photons coming towards 
us from the surface of last scattering. The peculiar velocity that we observe today, would 
be the sum of all string impulses with some suitable growth factors while the anisotropy 
would be due to the sum of the string impulses on the photons from the surface of last 
scattering. Therefore both processes depend on evaluating the magnitude of each impulse, 
the number of such impulses and, in the case of peculiar velocities, the growth of the 
velocity following the impulse. Since the calculations are so similar, the dependences on 
the string network parameters are also very similar. 

Let us first apply the multiple impulse approximation (Vachaspati T. 1992a; Perivolaropou- 
los L. 1993a,b) to obtain the mean and the standard deviation of cosmic string induced 
peculiar velocities. According to this approximation, the time between t eq (the time of 
equal matter and radiation) and the present (t p ) is divided into a set of Hubble time-steps 
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ti such that ti + i = 2tj. At every Hubble-step each long string gives a velocity impulse 
to matter within a distance £tj from it, where £ij is the radius of curvature of the string, 
and, by causality £ < 1. It is straightforward to count the total number of string induced 
impulses for each comoving scale Lh^Mpc. By including the growth factor (J^-) 1 /^ for 
the 11 th impulse, the velocity v{t p , L) of the scale L at t p may be written as: 

N L n s (L) 

v(t p ,L) = -AvJ2 knm(h 1/3 (1) 

n=0 m=l 

where 

Av = 57tG7^ s 7 s / 5 / 4 (2) 

and the indices n, m count Hubble steps and long strings within a Hubble volume respec- 
tively. In (1), (2) and in what follows we use the following notation: N L is the number of 
Hubble time steps during which a volume of comoving size L experiences coherent string 
impulses, n s (L) is the number of string impulses that the volume experiences per Hubble 
time, v s is the string velocity and 7 S the corresponding Lorentz factor, (£^0)6 is the bare 
string tension in units of 1CP 6 , 
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and T, fx are the renormalized string tension and string mass density. Simulations (Allen 
B. & Shellard E. P. S. 1990; Bennett D. P. & Bouchet F. R. 1988) indicate T ~ 0.7^ and 
[i ~ 1.4^0 where (xq is the bare string tension. The expression for / 5 ^ 4 in eq. (3) can be 
thought of in the following way: if we multiply / 5 ^ 4 by Av, the first term is the induced 
velocity due to the conical deficit of the string metric while the second term (proportional 
to 1/v) is due to the Newtonian potential of the wiggly string. In what follows, we shall 
not evaluate / 5 m using the simulation results but treat it like a free parameter. 

If we imagine a volume of comoving size L, the whole volume will get a coherent 
impulse due to strings moving outside the volume but within a distance equal to the 
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curvature radius of the string. If a string passes through the volume, it will not cause a 
coherent motion of the volume and so we neglect the effect of such strings. On the other 
hand, a wiggly string which is distant from the volume would still have some Newtonian 
gravitational effect for £ < 1 (expressed by second term in ^5/4) on the volume and, strictly 
speaking, we should take the distant strings into account with an effectively smaller value 
of ^5/4. (We thank Alex Vilenkin for emphasizing this point.) However, since we are 
keeping the number of strings (n s ) and / 5 / 4 as free parameters, the exact values do not 
concern us and we should view / 5 / 4 as an average value that accounts for all the strings. 
In addition, for £ ~ 1, the curvature radius is comparable to the horizon and our method 
automatically accounts for all strings. The number of strings that give coherent impulses 
to the volume can now be found: it is simply the number of strings within the horizon 
multiplied by the probability that the string passes within a distance ^t n from the volume 
of size L. Therefore, 

'& n — L(t r 



n s (L) = n s [ j (4) 

provided L(t n ) < ^t n and n s (L) = otherwise. Here L(t n ) = L(t n /t p ) 2 ^ is the physical 
size of the volume at the n th Hubble step and n s is the number of strings within the 
horizon. 

Now we square (1) and take averages with the assumption that the k vectors are 
random. That is, 

< k nm ■ k n > m > >= S nn iS mm i (5) 

This gives (Vachaspati T. 1992a) the present average (rms) peculiar velocity magnitude 
on a comoving scale Lh~ 1 Mpc: 

v(t p , L > it eq ) = 44a£/ 5 / 4 LJ 1 km/s (6) 

where 

« = \fn^{G^o)Q{v s ^ s ) , (7) 
5 



LlOO = L / (lOOh -1 Mpc) and we have taken » 1. The result (6) is valid on scales 
Lh~ l larger than the curvature radius of strings at t eq = 10h~ 2 Mpc because in deriving 
it we assumed L(to) > it eq . For scales L < £t eq we get 

L 



v(t p , L < Zteq) = 437 



1 + 2.9 1 - 



1/2 

af^i^h km/s . (8) 



The dispersion a(L) of v(L) may also be obtained (Vachaspati T. 1992a) as: 

a(t p ,L) = 0.36 v(t p ,L) . (9) 

As discussed above, the free parameter a involved in the peculiar velocity calculation 
can be pinned down by demanding that the cosmic string model be consistent with the 
recent detection of MBR anisotropy by COBE. The multiple impulse approximation can be 
used (Perivolaropoulos L. 1993a) to express the string induced MBR anisotropies in terms 
of a. In this case the approximation is realized as follows: A given photon experiences 
the impulses due to the various strings along its way from the last scattering surface to 
us. Each impulse can either red shift or blue shift the photon, depending on the various 
possible directions of the string velocity, photon momentum and string orientation. The 
temperature variations at some fixed angular separation will be correlated only when the 
photon beams at that angular separation will experience the same impulse. The correlation 
function is thus obtained by counting the number of common impulses on each angular 
scale 9. In this case, due to the lack of a growth factor, the dependence on £ (the string 
curvature radius) comes only as a square root dependence in counting the number of strings 
that come within a radius of curvature distance from the photon. From the correlation 
function it is straightforward to obtain the rms temperature fluctuation for angular beam 
separation equal to 60° and smoothing on a scale of 10° as in the COBE experiment. The 
result is (Perivolaropoulos L. 1993a): 

(^) =1.6axl0- 5 (10) 
\ 1 / rms 
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where a is the same product of parameters as in (7). 

Now the COBE observations fix the rms temperature fluctuations to be (^j^)rms = 
(1.1 ± 0.2) x 1CP 5 which when inserted in (6) gives 

a -0.7 ±0.2 (11) 

It is of interest to use the values of evolution parameters n s , v s ^ s and £ obtained in 
numerical simulations to extract from the value of a in (11), a rough value of the single 
free parameter of the string model G/j,q. From (Allen B. & Shellard E. P. S. 1990) we have 
n s — 10, v s — 0.15 and £ — 0.7. Using these values in (7) we obtain (Perivolaropoulos L. 
1993a): 

G/j = (1.8 ±0.5) x 10~ 6 (12) 

which is the value required from galaxy and large scale structure calculations (Perivolaropou- 
los L., Brandenberger R. & Stebbins A. 1990; Turok N. & Brandenberger R. 1986; Sato 
H. 1986; Stebbins A. 1986). (The value of G\i given in (Perivolaropoulos L. 1993a) refers 
to the effective Gfi which includes the string small scale structure and is therefore slightly 
different from (12)). The same value for G\i has been obtained in (Bennett D. P., Stebbins 
A. & Bouchet F. R. 1992) by using the numerical simulations of (Bennett D. P. & Bouchet 
F. R.1988; Bouchet F. R., Bennett D. P. & Stebbins A. 1988) to simulate the MBR sky 
and compare with the COBE results. This is an additional test of consistency for our 
results. 

Now, using (11) in (6) and taking into account a la spread (eq. (9)) we get 

V(t p , L > it eq ) = (31 ± 11) C/s/^ioO km / S ( 13 ) 



v(tp,L < £teq) = (306 ± 110) 



1 + 2.9^1-— J j f 5/4 hkm/s. (14) 



The dependence of v on the length scale is plotted in Fig. 1 for / 5 ^ 4 = 1.8, £ = 1 and 
h = 0.5. The observations (Burstein D. et. al. 1983; Collins C, Joseph R. & Robertson 

7 



N. 1986; Dressier A. et. al. 1987; Groth J.,Juszkiewicz R. & Ostriker J. 1989) are also 
indicated with their error bars. 

The figure shows that the peculiar velocities from cosmic strings are in good agreement 
with observations on small scales but not on large scales. The inconsistency is not an 
artifact of our choice of parameters. For example, we could take £ = 1 (the largest 
permissible value), h = 0.5 (the smallest value of h) and / 5 / 4 = 7 so that the observation 
on 60/i _1 Mpc is barely accommodated within the la prediction of the string model. But 
then, the predicted la velocity on scales of 8/i _1 Mpc range from 1135km/s to 2409km/s. 
These velocities are clearly inconsistent with observations. In other words, no choice of 
parameters will allow us to explain the velocity observations on all scales. 

We now consider the Cosmic Mach Number (CMN) induced due to strings. The CMN 
has been shown to be independent of various uncertainties such as biasing and hence 
may be a valuable statistic in comparing theory with observation (Ostriker J. & Suto Y. 
1990). As we shall see, the string induced CMN only depends on the parameter £ and for 
reasonable values of £ the agreement between the string predictions and observations on 
small scales is quite good. Hence we believe that cosmic strings can successfully explain 
the observations on small scales but fail in explaining the large scale data. 

The Cosmic Mach Number M(L, a) is defined (Ostriker J. & Suto Y. 1990) as the ratio 
of the rms velocity v(t p ,L) of a volume of size Lh^ 1 over the dispersion of the velocity 
field smoothed on a scale a, in the rest frame and within the same volume: 




v(t p ,L) 



(15) 




Inserting (8) for the peculiar velocities gives us (for a, L < £,t eq ), 




(16) 
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If a < ^t eq < L, we get, 

/ 7£ 2 t 2 

M ' L ' a >=V 2.7^-2^-0.7^ - (17> 

Observations on scales larger than 5/i _1 Mpc indicate (Ostriker J. & Suto Y. 1990) that 

M(L = Sh~ l Mpc, a = hh^Mpc) = 2.2 ± 0.5 

(18) 

M(L = 18/i _1 Mpc, a = 5h~ 1 Mpc) = 1.3 ± 0.4 . 
while the corresponding predictions of the cosmic string scenario obtained from (16) for 
£ = 1, are 

M S (L = 8/i _1 Mpc, a = hh~ l Mpc) = 2.5 

(19) 

M S (L = ld>h~ l Mpc, a = 5/i _1 Mpc) = 0.8 . 
The CMN predictions are smaller for smaller values of £. For example, if £ = 0.7, one gets 

1.9 and 0.5 in eq. (19). Simulations indicate £ ~ 0.7, and so we can say that the cosmic 

string predictions are in rough agreement with the CMN observations. If, however, the 

simulations can pin down £ better, the CMN observations would be a powerful test for the 

cosmic string scenario. 

The success of the string scenario with the observations on small scales suggests that 
we should choose our parameters so that both the peculiar velocity magnitudes and the 
CMN are in good agreement on these scales. That is, we choose £ = 1> /s/4 = 1-8 and 
h = 0.5. With this choice of parameters, the peculiar velocity on 60/i _1 Mpc is 89±30km/s 
and is off from the observations of roughly 600km/s by a factor of about 7. This gives us 
an idea of the magnitude of the problem. 

We are unable to give a spread on the predicted values of the CMN because our analytic 
model does not take the correlations of velocity impulses on the large and small scales into 
account. These correlations are important for the CMN because the dispersion occurring 
in the denominator of eq. (15) is the dispersion of the smaller scale volumes (size a) within 
a bigger volume (size L). For example, it is quite possible that if the bigger volume has a 
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larger than average velocity, then the smaller volumes will also have a larger than average 
velocity. Such correlations cannot be taken into account in our model since we can only find 
the average velocities of the small and big volumes. It would also be wrong to include the 
spread in peculiar velocities (eq. (9)) and simply propagate these into the CMN because 
the calculated spread includes all volumes of size a and not just those contained within 
some larger volume of size L. 

The prediction of the cosmic string model is in better agreement with observations 
than adiabatic CDM which predicts (Ostriker J. & Suto Y. 1990) M(L, a = hh~ x Mpc) to 
be well below 1 on both 8h~^Mpc and 18h~^Mpc. 

In conclusion, we have used the multiple impulse approximation to identify two new 
features of the cosmic string model valid with both hot and cold dark matter: 

a) The Cosmic Mach Number predicted by the model is consistent with observations 
on scales 5 — 20/i~ 1 Mpc. 

b) The magnitude of the predicted large scale velocity flows normalized by the COBE 
detection of MBR anisotropy is consistent with observations on small scales for reasonable 
values of the string evolution parameters but the large scale peculiar velocities cannot be 
explained with these parameters. 

The potential problem pointed out in b) is also encountered in other theories like 
standard CDM where a similar L _1 scaling of the magnitude of large scale velocity flows 
has been shown (Kaiser N. 1983; Vittorio N. & Silk J. 1985; Vittorio N. & Turner M. S. 
1987; Kolb E. W. & Turner M. S. 1990). 
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FIGURE CAPTIONS 



1. The predicted peculiar velocity versus the scale lies in the shaded area between the 
two shown curves within a la error including the standard deviation of v. The data 
points with error bars are the observations. The parameter values / 5 / 4 = 1.8, £ = 1, 
h = 0.5 were used. 
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